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Abstract 

We argue that topological Kondo insulators can also have ‘intrinsic’ topological order associated with 
fractionalized excitations on their surfaces. The hydridization between the local moments and conduction 
electrons can weaken near the surface, and this enables the local moments to form spin liquids. This 
co-exists with the conduction electron surface states, realizing a surface fractionalized Fermi liquid. We 
present mean-field solutions of a Kondo-Heisenberg model which display such surfaces. 
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I. INTRODUCTION 


An important development of the past decade has been the prediction and discovery of topologi¬ 
cal insulators (TI) [1-7]. These materials are well-described by traditional band theory, but possess 
strong spin-orbit interactions that result in a non-trivial winding of the ground state wavefunction 
in a manner analogous to the integer quantum Hall effect. Since their discovery, the multitudinous 
effects of interactions have been a prominent topic of study. One compelling proposal to emerge 
is the notion of a topological Kondo insulator (TKI) [8-10]. In contrast to a band insulator, a 
Kondo insulator only develops an insulating gap at low temperatures, and the magnitude of the 
gap is controlled by electron-electron interactions. Doniach explained this phenomenon through 
the Kondo lattice model [11] in which a lattice of localized moments is immersed within a sea of 
conduction electrons. At high temperatures, RKKY-type exchange interactions dominate and an 
ordered magnetic state results. Conversely, at low energies, strong interactions between localized 
moments and conduction electrons become important; the system crosses over into either a metal¬ 
lic phase well-described by Fermi liquid theory (FL) or, if the chemical potential is appropriately 
tuned, a Kondo insulator. As strong spin-orbit coupling is often present in these materials, the 
possibility that a Kondo insulator may have a nontrivial topological character is well-justihed. 

Of specihc interest has been the Kondo insulator samarium hexaboride (SmBg). A number of 
experiments have examined the proposal that it is a TKI; transport measurements have established 
the presence of metallic surface states [12-16], and angle-resolved photoemission spectroscopy 
(ARPES) results appear consistent with the expected Dirac surface cones [17-21]. Nonetheless, 
the spin-polarized ARPES measurements [21] remain controversial. 

However, as the TKI phase is well-described within a mean held framework [10], its topological 
properties are not expected to be markedly different from what has already been observed in its 
uncorrelated cousins. More intriguing is the potential the topologically protected surface states 
present for new interesting phases [22-25]. In SmBg, this expectation is motivated experimentally 
by ARPES measurements which End light surface quasiparticles [17-19] in contradiction to current 
theories which predict heavy particles at the surface [8, 9, 26, 27]. Ref. [24] proposes “Kondo 
breakdown” at the surface as an explanation. They show that the reduced coordination number 
of the localized moments at the surface may lead to a suppressed Kondo temperature. At low 
temperature these moments are thermally decoupled from the bulk. 
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In this paper, we propose the existence of a fractionalized Fermi liqnid (SFL*) on the surface of a 
TKI. This state is characterized by ^''intrinsic topological order” on the surface of a TKI, in which 
the local moments form a spin liquid state which has ‘fractionalized’ excitations with quantum 
numbers which cannot be obtained by combining those of one or more electrons [28] . Rather than 
being thermally liberated, as in Ref. 24, the surface local moments exploit their mutual exchange 
interactions to decouple from the conduction electrons, and form a spin liquid state, as in the 
fractionalized Fermi liquid state (FL*) [29, 30]. We will present mean-held computations on a 
Kondo-Heisenberg lattice model which demonstrate the formation of mutual singlets between the 
surface local moments, while conducting surface states of light electronic quasiparticles are also 
present. 

Somewhat confusingly, our SFL* state is ‘topological’ in two senses of the word, a consequence 
of unfortunate choices (from our perspective) in the conventional terminology. As in conventional 
TI, it is ‘topological’ because it has gapless electronic states on the surface induced by the nature 
of the bulk band structure. However, it is also ‘topological’ in the sense of spin liquids [28], because 
of the presence of fractionalized excitations among the local moments on the surface. 

The outline of our paper is as follows. We specify our Kondo-Heisenberg model in Section H. 
In Section HI, we present the mean-held solution of this model for the case of a translationally- 
invariant square lattice with periodic boundary conditions. The ehect of the surface on the mean 
held solutions is addressed in Section IV where the presence of the SFL* state is numerically 
demonstrated. We conclude in Section V with a discussion of our results and their relevance to 
physical systems. 


II. MODEL 

Here we present the specihc form of the Kondo-Heisenberg lattice model to be studied: 

H = H, + Hh + Hk. ( 1 ) 

The hrst terms represents the hopping Hamiltonian of the conduction electrons, 

Hc = -Y^ tij + h.c.) , (2) 

(o> 
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where the operator creates an electron at site rj of spin a =■!■, -I- The remaining two terms 
establish the form of the interactions: Hh is a generalized Heisenberg term which specihes the 
inter-spin interaction while Hk is a Kondo term and describes the electron-spin exchange. 

The spin-orbit conpling of the /-orbital imposes a classihcation in terms of a (2 J -|-1) mnltiplet, 
where J is the total angnlar momentnm. In general, this degeneracy is farther lifted by crystal 
helds and we will consider the simplest case of a Kramers degenerate pair of states. We start from 
an Anderson lattice model [31] with hopping tf between /-orbitals and onsite repulsion Uf. To 
access the Kondo limit {Uf —)■ cxd), we perform a Schrieffer-Wolff transformation [32] and obtain a 
term of the form 


(o> 


( 3 ) 


where fia creates a spinon at site rj, and Jh ~ t) /Uf. This limit imposes the constraint 
Yla fiafioi = ^ further ensures that the correct commutation relations for the “spin” oper¬ 
ators S/ = ci-re obeyed. By using the Fierz identity (and dropping a constant) we can 

verify that we indeed have the familiar Heisenberg term: 




1/3/2/? ■ ff js Jh ^ ^ Sj ■ Sj , 


( 4 ) 


(o> 


(o> 


where cr = (a^, It is important to note that the spinon operators fja do not have a 

uniquely dehned phase. In fact, by choosing to represent the spins in terms of constrained fermion 
operators, we are formulating the Kondo lattice model as a U(l) gauge theory. This emergent 
gauge structure is what permits a realization of the fractionalized phases we will discuss [29, 30]. 

For the electron-spin interaction, Hk, we follow the construction of Coqblin-Schrieffer [33] for 
systems with spin-orbit coupling. In order for the interaction to transform as a singlet, the electron 
and spin must couple in a higher angular momentum channel. For simplicity, we assume a square 
lattice and that the spins and conduction electrons carry total angular momentum differing by 
/ = 1. In the Anderson model, an appropriate interaction term is then 


Hint ~ (asin -i sin ky) \ f) {f^;a\ h.c. (5) 

k,a 
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For instance, the interaction between moments with total angular momentum J = 3/2 and spin- 
1/2 electrons would take this form. We will verify in the next section that for the purpose of 
obtaining a TKI, this coupling is sufficient. We next dehne the electron operators 

dka 2 (cr sin fca; -|- isinky') Cko, dia ^®(Cj+x,Q Ci—x^a) (Ci+y,a Cj—y,a) • (6) 

and, taking the same Uf ^ oo limit as above, again implement the Schrieffer-Wolff transformation 
[32] to obtain 

i 

where Jk ~ V^/t//. 


We next perform a Hubbard-Stratonovich transformation of the Kondo and Heisenberg terms: 


H' =H, + Ho 

+ h.C.) + ^ X 

(ij) 



X! \xjfJ'fj+p.,afjo^ + Xjy.fjafj+p.,a + ^jfjafja 


ja/j. 
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\^ 3 ? , IwfI' 
Jk - ‘^Jh 

IJ- 


( 8 ) 


We proceed with a saddle-point approximation, and treat the helds Vj, and Xj as real constants 
subject to the self-consistency conditions 




1 = 



XjfJ. Jh (^fjafj+P-,^^ y 


(9) 

( 10 ) 


This can be formally justihed within a large-iV expansion of Eq. 1, with N the number of spinons. 
As we are specifying to the case of an insulator, it further makes sense to require perfect half-hlling. 
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Since Uf = 1 already, this results in a final equation for the chemical potential 
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(11) 


III. TRANSLATIONALLY INVARIANT SYSTEM 


We begin by solving Eqs. 9 — 11 in a translationally invariant system with periodic boundary 
conditions on a square lattice. Letting Vj = V, Xjx = Xjy = X) and = p, we perform a 

Fourier transform: 


k 


H(k) 


fhik) 0 \ 
\ 0 h*i-k)) 



/kf /kj.) 


h(k) = 


ec(k) V{sin kx + i sin ky) 

V{sin k^ — i sin ky) ^/(k) 


( 12 ) 

(13) 


For simplicity, we only consider nearest-neighbour coupling between spins; for the electron dis¬ 
persion, a slightly more general description is required and we also take next-nearest neighbour 
hopping into account. The dispersions are given by 


ec(k) = —fi(cos kx -l- cos ky) — 2^2 cos kx cos ky — /i, e/(k) = —x(cos kx + cos ky) + A (14) 

where the subscripts “c” and “/” refer to the electrons and spinons respectively. In the following, 
we will use units of energy where ti = 1.0. 

Since TI’s exist as a result of a band inversion, it’s important to ask which sign y; will take. 
Naturally, when V = 0, the particle-hole symmetry of our mean field ansatz implies that x > 0 
and X < 0 have the same energy. At finite hybridization, however, one will become preferable. We 
note that when x and ti have opposite signs, the energy of the lower band will be less than the 
Fermi energy and hence occupied throughout most of the Brillouin zone (BZ): an increase in V 
will push most of these states to lower energies. Conversely, if x and ti take the same sign, in one 
of part of the BZ no states will lie below the Fermi energy while in another both the upper and 
lower band will. It therefore makes sense to expect sign (x) = —sign {ti). In the parameter regime 
explored, the numerics always find this to be the case. 

By construction, the Hamiltonian Hi supports a non-trivial topological phase and is in fact the 
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familiar Bernevig-Zhang-Hughes model [4] used to describe the quantum spin Hall effect in HgTe 
wells. We can see this by studying the eigenfunctions of h(k); 


V’±(k) 


1 

'\/2d(cfi~d3y 


^3 i d 

V (sin kx + i sin ky) 


( 15 ) 


where d(k) = \/+ l/^(sin^ k^ + sin^ ky), and d 3 (k) = (ec(k) — e/(k)) /2. If d 3 (k) > 0 or 
d 3 (k) < 0 for all k, these functions are well-dehned on the entire BZ and the system is in a 
topologically trivial phase [34], If this is not the case then it is impossible to choose a globally de- 
hned phase - the ground state wavefunction has nontrivial winding and characterizes a topological 
insulator. From Eq. 14, we see that this occurs when 


-2 < 


/i + A + 2^2 
ti-X 


< 2 . 


(16) 


Alternatively, we can obtain the same result by calculating the Z 2 invariant u [35]: when Eq. 16 
holds, u = —1 and the system is a TI. 


We will typically be studying systems with |ti| 3> |x| and 1^2/Hi small (implying fi and A small 
as well), so Eq. 16 is not difficult to fulhll. In Fig. 1 the energy spectrum of the system in a 
slab geometry is shown for Jh = 0.15, Jk = 0.3. Half-hlling is maintained on every site (see 
Appendix A), but V and y were determined by self-consistently solving Eq. 9 in a periodic system. 
In Fig. 1(b), the topologically protected Dirac cone is clearly visible. 


If we ignore the effect the boundary will have on the values of V, y, and A, we can calculate 
the Fermi velocity of the Dirac cone [34]: 


vp — 2V 



1 ~ 2^2)1 
2H| + |y| 


~ V\x\ 


(17) 


where we’ve assumed |y| <C in the second equation. This is consistent with the prediction that 
the quasiparticles at the surface be heavy [8, 9, 26, 27]. For the parameters shown in Fig. 1, this 
formula predicts vp = 0.0592, consistent with the numerically determined value vp = 0.0585. 
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FIG. 1. Energy spectrum for Jh = 0.15, Jk = 0.3. Both V and y are constant throughout the bulk, 
but both /Ji and have been self-consistently solved to ensure that Uc = n/ = 1 on every site (see 
Appendix A), (a) The full spectrum is shown, (b) A closer view of the insulating gap, where the Dirac 
cone is clearly visible. We use units with ti = 1.0. Calculations were done with t 2 = —0.25 and at a 
temperature of 10“®. 

IV. SYSTEM WITH BOUNDARY 


We now consider the effect the boundary will have on the mean held conhguration and demon¬ 
strate the presence of two new fractionalized phases. Generally, we expect that the lower coordi¬ 
nation number at the boundary will suppress the (nonlocal) hybridization: 14urf ~ 3Vbuik/4. While 
the decrease in 14urf will induce an increase in the spinon bond parameter Xifj. [36] both parallel 
and perpendicular to the surface, the parameter parallel to the surface will be more strongly af¬ 
fected. Since Heisenberg coupling ultimately favours an alternating bond order, in the absence of 
hybridization V, this anisotropy will result in a further decrease in the magnitude of the spinon 
bond parameter perpendicular to the surface, |x_l|- 

When these effects are predominant, an FL* on the surface is realized: the hybridization V) 
vanishes on one or more layers at the surface and x± vanishes on the innermost layer. The existence 
of the SFL* phase is shown numerically by self-consistently solving Eqs. 9—11 in a slab geometry 
and comparing ground states energies (some details are given in Appendix A). The resulting phase 
diagram is shown in Fig. 2(a). In fact, we hnd two distinct SFL* phases: a decoupled spin chain 
and a decoupled spin ladder, which are depicted in Figs. 2(b) and 2(c) respectively. In Fig. 3 we 
plot the spatial dependence of the mean held parameters in both SFL* states. The plots in the 
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(a) 


(b) 


(c) 


Kondo coupling 

FIG. 2. (a) Schematic phase diagram of surface states. (b),(c) Cartoon depictions of surface FL* states. 
In the dark blue region, the electron spins and localized moments are locked into singlets. Towards the 
edge (the pale blue region outlined in orange) the conduction electrons decouple from the moments, and 
the latter form a spin liquid. Naturally, the conduction electrons remain coupled to each other at all sites. 

left column correspond to a spin chain SFL* state whereas the right column corresponds to a spin 
ladder SLF* state. The phases are distinguished by whether V vanishes on the hrst site only or 
on both the hrst and second site, shown in Fig. 3(a) and (b) respectively. In Figs. 3(c) and (d) 
our intuition regarding the behaviour of y near the boundary is conhrmed: |x_l| is suppressed to 
zero whereas |x||| increases to the value it would assume in a single dimension. The huctuations 
of the Lagrange multiplier held A (Figs. 3(e) and (f)) are a rehection of the on-site requirement of 
half-hlling for both the spinons and electons. 

In Fig. 4(a), the spectrum of the spin chain SFL* state is shown. The red dash-dotted curve is 
the dispersion of the spinons calculated at mean held. While we do not claim that this accurately 
represents the Heisenberg chain, we nonetheless expect gapless spin excitations [37]. The remaining 
in-gap states can be understood as the result of the mixing of the surface layer of conduction 
electron with the Dirac cone. Consistent with its topology, even if the Dirac cone is no longer 
present at the chemical potential, two chiral bands traverse the gap from the conduction to the 







\' 
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V 
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decoupled spin ladder 
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FIG. 3. Spatial dependence of mean field parameters in SFL* phases. In the left column, we plot values 
corresponding to the spin chain SFL* {Jh = 0.15, Jk = 0.3) while on the right values corresponding to 
the spin ladder SFL* {Jh = 0.25, Jk = 0.3) are shown. (a),(b) Hybridization V). (c),(d) Spinon bond 
parameters Xifj, the direction perpendicular (blue) and parallel (red) to the boundary. (e),(f) The La¬ 
grange multiplier field A*. In (a) —(f), the yellow dashed line plots the value obtained in the translationally 
invariant case. We use units with ti = 1.0. Calculations were done with t 2 = —0.25 and at a temperature 
of 10“^ 10 


































FIG. 4. Energy spectra in SFL* phases, (a) Spin chain SFL* {Jh = 0.15, Jk = 0.3). The ground state 
has V) = 0 on the hrst surface layer and the moments form a spin chain decoupled from the bulk, (b) 
Spin ladder SFL* {Jh = 0.25, Jk = 0.3). The ground state has V) = 0 on the first two layers and a spin 
ladder is present on the surface. In both figures, the dash-dotted red curve represents the one-dimensional 
cosine dispersion found for the spinons and is merely an artifact of the ansatz. We use units with ti = 1.0. 
Calculations were done with t 2 = —0.25 and at a temperature of 10“^. 


valence band and the snrface is metallic. In this case, an additional four metallic surface states 
per spin are present, but these are not topologically protected and we can imagine pushing them 
below the chemical potential in a number of ways, such as, for instance, softening the restriction 
imposed by Eq. 11. 

The spectrum corresponding to the second surface FL* state, the decoupled spin ladder, is 
shown in Fig. 4(b). The red curve representing the spinons is now two-fold degenerate per spin (a 
small splitting is hidden by the thickness of the line). Even more so than for the spin chain, this 
result is an artifact of the mean held calculations: save in the limit where the legs of the ladder 
are completely decoupled, a ladder of spin-1/2 particles is gapped [37, 38]. 

In both phases, the metallic bands have lighter quasiparticles than predicted by the transla- 
tionally invariant theory in Eq. 17. For the spin chain, the surface velocity of the leftmost state 
in Fig. 4(a) is vp = 0.095, compared to vp = 0.052 for the Dirac cone of Fig. 1(b). For the 
spin ladder, the effect is even more pronounced. There, the lightest state has a Fermi velocity of 
vp = 0.48 compared to the translationally-invariant value vp = 0.072. 

The structure of the fractionalized excitations in the SFL* states found here is rather simple: 
just a free gas of neutral 5 = 1/2 spinon excitations. We view this mainly as a ‘proof of principle’ 
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that such SFL* states can exist on the surface of TKI. Clearly, more complex types of spin liquid 
states are possible on the surface, and also in three-dimensional TKI with two-dimensional surfaces. 


V. DISCUSSION 


The strong electron-electron interactions in topological Kondo insulators make them appealing 
candidates for searching for novel correlated electron states. In many heavy-fermion compounds, 
the strong interactions acting on the /-electron local moments are quenched by the Kondo screening 
of the conductions electrons, and the resulting state is eventually a Fermi liquid, or a band insulator 
for suitable density. The topological Kondo insulators offer the attractive possibility that the 
hybridization between the local moments and the conduction electron states can be weakened near 
the surface [22, 24], and this could explain the light effective masses associated with the surface 
electronic states [17-19]. With the weakened hybridization, we have proposed here that the local 
moments may form a spin liquid state with ‘intrinsic’ topological order. As the fractionalized 
excitations of such a spin liquid co-exist with the conduction electron surface states similar to 
those of a conventional TI, the surface realizes a fractionalized Fermi liquid [29, 30]. 

This paper has presented mean-field solutions of Kondo-Heisenberg model on a square lattice 
which act as a proof-of-principle of the enhanced stability of the such surface fractionalized Fermi 
liquids (SFL*). We fully expect that such solutions also exist on the surfaces of three-dimensional 
lattices, relevant to a Kondo insulator like SmBg. 

The recent evidence for bulk quantum oscillations in insulating SmBg [39] is exciting evidence 
for the non-trivial many-electron nature of these materials. It has been proposed [40] that these 
oscillations appear because the magnetic-held weakens the hybridization between the conduction 
electrons and the local moments, and this releases the conduction electrons to form Fermi surfaces 
leading to the quantum oscillations. The fate of the local moments was not discussed in Ref. 40, 
but a natural possibility is that they form a bulk spin liquid, similar to the surface spin liquid we 
have discussed here. Thus, while we have proposed here the formation of a SLF* states in SmBg 
in zero magnetic held, it may well be that a bulk FL* state forms in high magnetic held. 
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Appendix A: Mean field theory with bonndary 


In this appendix, we consider the mean held eqnations in the presence of a boundary. We dehne 
the lattice to be hnite in the x-direction, Xi = 1,..., N, and inhnite in the ^/-direction (to remove 
factors of i, we actually switch the x- and ?/-directions compared to Eq. 6). We rewrite the Fourier 
transform, which is now only valid in the ^/-direction: 



Translational invariance in the ^/-direction implies that the mean held parameters will depend 
only on the distance from the boundary - the roman indices i, j, etc. label the x-coordinate only. 
We express the Hamiltonian in block form as 


k 


n,{k) = 




t ^ 

k<7 


■t t 

l/ccr’ ^2ka’> * 




t 

1 /C(T’ 




2ka'! 


= 




l+N ^k(7'> ^2+N 




t 


H^{k) = H^{-ky 


hc{k) hcf{k) 

hi,(k} h,(k) 
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(A2) 


with blocks given by 



2Vism.k —V 2 0 

Vi 2V2sin/c — V3 
0 V 2 2 V 3 sin k 

To determine correlation functions, we diagonalize the Hamiltonian numerically. For each k, we 
hnd the matrices U{k) such that 

U^n{k)U{k) = A{k), Aij{k) = 6ijEj{k). (A 6 ) 
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Then, the mean held equations of Eqs. 9 — 11 may be expressed as 


(A7) 


i=l 


i = E 


dk 


2N 

= E / § = 2 / ^ E n(E,(k))U„(k)Ul(k) (A8) 

a ^ 1=1 


J2 [ ^ {4kaUa) = / T“ f^asin A; /4„/ifca) + 


J, 




dfc 


27r 

dk 

2Ar 


/ TT hasinfc /V’L^i+iV,fea) + (V’l+pto^i+iV.fea) - (V’Ll,fcaV'i+iV,fca 


i=l 


-Jk I — Y,<Ei{k)) 2sinkU,+NAk)UlM + U^+NAk)Ul^^{k) -U,+NAk)Ul_^{k) 


(A9) 


Jj^Y^ 

Aix ^ ^ 


2 

Jh 


E 


dk 


‘ dk r 

‘ dk 

^ L 

2N 


fi+l,kaf{fikafk+^t 


,ka 


'^i+N+l,ka'^i+N,ka ,ka 


■’« I ^E”(^'('=)) C'i+w(*^)f7;.+„+i(*:) + A+A.aj(*:)di+«+i('=) 


i=l 


Xij/ — 


Jh 


5^ / ^2 cos k i^flkJika) = JhY^ / ^ ^ (^/'!+iV,fcaV'i+iV,fca 


dk 


2N 


2* / l^T.'^(E,(k))m!ikU,+^,,(k)U;,^M 


1=1 


(AlO) 
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